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Abstract 

We derive a partially gauge fixed Hamiltonian for black hole formation via 
real scalar field collapse. The class of models considered includes many the- 
ories of physical interest, including spherically symmetric black holes in D 
spacetime dimensions. The boundary and gauge fixing conditions are cho- 
sen to be consistent with generalized Painleve-Gullstrand coordinates, in 
which the metric is regular across the black hole future horizon. The result- 
ing Hamiltonian is remarkably simple and we argue that it provides a good 
starting point for studying the quantum dynamics of black hole formation. 
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1 Introduction 



It is widely believed that black hole thermodynamics provides important 
clues about the structure of the underlying quantum theory of gravity. Spher- 
ically symmetric black holes provide a rich testing grounds for analyzing these 
clues. Moreover, recent work pQ El E] suggests that despite the absence of 
gravitational dynamics, the spherically symmetric sector may contain suffi- 
cient structure to describe the origin of black hole entropy. 

Perhaps the most significant feature of black hole thermodynamics is 
its apparent universality. The famous Bekenstein formula for black hole 
entropy seems to be valid for a very large class of black holes, irrespective 
of the underlying classical or microscopic quantum dynamics. It is therefore 
important to understand the classical and quantum structure of black holes 
generically, and to find new approaches to analyze them. 

The purpose of this letter is to present the results of a novel Hamiltonian 
analysis that may be very useful in the analysis of spherically symmetric black 
hole formation via scalar field collapse. The general formalism we use is based 
on 2-D dilaton gravity but the class of models we consider describes many 
specific black holes of physical interest, including the spherically symmetric 
sector of higher dimensional Einstein gravity. 

In order to have a system of equations that is valid at, or even below, 
any event horizons that might form, we look to Painleve-Gullstrand (PG) 
coordinates. These were used several years ago by Kraus and Wilczek [I] to 
derive a self-interaction correction to black hole radiation in the collapse of 
spherical shells. Subsequently, Louko et al[5\ analysed in detail the Hamil- 
tonian dynamics of spherical shell and dust collapse using PG boundary 
conditions. More recently PG coordinates were used in Hamiltonian analysis 
of spherically symmetric collapse in 4-d Einstein gravity [6], and to examine 
the quantization of generic single horizon black holes[71 [8]. 

A crucial ingredient in our analysis is that we do not completely fix the 
gauge. As in [3[8], only the spatial diffeomorphisms are fixed, so that one is 
left with a midi-superspace model containing a single Hamiltonian constraint 
that generates time diffeomorphisms. The proposed analysis of the quan- 
tum dynamics will therefore require imposing this Hamiltonian constraint 
on physical states and then solving for the time evolution generated by the 
remaining Hamiltonian. The latter is non-vanishing even on the constraint 
surface because of the presence of the usual boundary terms needed to make 
the classical variational principle well defined. Although completely reduced 
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models are in principle simpler due to the absence of an infinity of unphysical 
modes, the resulting Hamiltonians tend to be very complicated, highly non- 
linear and non-local [9j. Moreover, since the degrees of freedom that account 
for black hole entropy may reside in diffeomorphisms that become physical 
in the presence of black hole boundary conditions [1] one may be "throw- 
ing out the baby with the bath water" by completely gauge fixing prior to 
quantizing. 

The most striking result of our derivation is the simplicity of the final, 
partially reduced Hamiltonian. It consists of three terms. The first is purely 
gravitational, being proportional to the mass observable. The second term 
is essentially the matter energy density, which contains only scalar field vari- 
ables. Finally there is a single term that couples the gravitational mode to 
the scalar degrees of freedom. The coupling is a fourth order polynomial in 
the phase space variables, and only second order in the momenta. We will 
argue that this Hamiltonian system provides an excellent starting point for 
studying the quantum dynamics of black hole formation. 

2 Classical Action and Solutions 

The formalism we will use to describe generic single horizon black holes is 
that of 2-D dilaton gravity (for a review see [10J). This approach has been 
successfully used to study spherically symmetric scalar field critical collapse 
in any dimension The action is: 

Sa\g, = 2^/ dxdty/=g (^R(g) + . (1) 

where I is a parameter with dimensions of length which is generally taken to 
be the Planck length in the theory. The matter action is: 

S M bP,g,<f>] = -\ J ' d 2 x^gh^)\Vij\ 2 (2) 

Different models are distinguished by two functions of the dilaton field; the 
dilaton potential V(0) and the coupling h(4>) between the matter field and the 
dilaton. This class of models describes a large number of physically relevant 
theories, including the spherically symmetric sector of higher dimensional 
Einstein gravity. In the case of four dimensions, for example, <fi ~ r 2 is the 
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area of the invariant 2-sphere at radius r, V ~ l/\/0 an d ^ ~ 0- For details 
of the Hamiltonian analysis of the vacuum sector, see for example|12j. 

The vacuum theory has no propagating modes and is exactly solvable. 
There exists a one parameter family of static solutions, which can be repre- 
sented in Schwarzschild like coordinates x = as: 

ds 2 = - 2GMl)dt 2 s + {j{4>) - 2GMl)- 1 dx 2 (3) 

where j(0) = f$ d<pV(<f)) and M is the mass. 

This metric is not asymptotically flat, but, for spherically symmetric grav- 
ity in d = n + 2 spacetime dimensions, it is related to the (r, t) components 
of the physical, asymptotically flat higher dimensional metric ds 2 hy by the 
conformal reparametrization: 

,2 ./,n,2 -t,s( /- 2GML ,„ . 2GMI , 2 \ 
^ 2 = j(0)c^ = j(0) f-(l - -7^)^' + (1 " -J^-y dr ) ( 4 ) 

where the radial coordinate is related to the dilaton by: 

dr = l^t- oc (5) 
3 VP) 

Note that this conformal reparametrization corresponds to a local field re- 
definition which converts the action in Eq. (TjQ) to its more standard form of 
spherically symmetric gravity involving the physical metric and containing a 
kinetic term for the dilaton field. This field redefinition is singular at <fi = 0, 
but this surface must in any case be excluded from the manifold in order to 
represent black hole spacetimes. In fact = will be treated as a boundary 
at which suitable boundary conditions must be imposed. 

We assume that j(<ft) is monotonic, so that there is at most one horizon, 
and without loss of generality that it vanishes at = 0. These conditions 
guarantee the existence of an apparent horizon in the classical vacuum solu- 
tion. 

It will be useful to consider a coordinate system in which the physical 
metric is in PG form, i.e. 



ds 2 = j(0) [-dt 2 + (dx + f-^dt) 2 )j (6) 
where the time coordinate is related to the Schwarzschild time by: 



2GMI dx 
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In this coordinate system the physical (asymptotically flat) spatial metric 
is constant and the entire metric is regular at the horizon. The slicing can 
therefore be extended from spatial infinity all the way to the singularity at 
= 0. It is therefore well suited for the study of horizon formation and 
evolution. 



3 Hamiltonian Analysis 

An analysis of the vacuum theory has been presented in [12] . Here we present 
the canonical analysis of the full theory (i.e. with matter), using the metric 
as parametrized in modified ADM form: 

ds 2 = e 2p (-a 2 dt 2 + (dx + Ndtf) . (8) 

We assume that the 2-D spacetime manifold M 2 = R x H, where E is a 
spacelike slice. Note that the exponential parameterization for the conformal 
mode is traditionally used because it manifestly restricts the conformal mode 
to be positive definite. 

As usual the lapse and shift functions a and N, respectively, are lagrange 
multipliers that give rise to the Hamiltonian and diffeomorphism constraints 
in the theory. The momenta of the remaining fields are given by the following: 

p = Np' + N' -aGU^ (9) 
= N(f)' -aGTl p (10) 

« = (ID 

where the prime denotes partial differentiation with respect to the spatial 
coordinate x. The Hamiltonian is: 

H = Jdx(aG + NF) + H B (12) 
where Q is the Hamiltonian constraint, 

~G~~G~ p ~ 2G~P~ 1 ^ ' 



with 

i ( n 2 

n 



2 \h((f)) 
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and T is the diffeomorphism constraint: 



T = P 'u p - u' p + + if/Hp (15) 

Hb is a boundary term required to make the variational principle well defined. 
It will be determined by the boundary conditions on the fields, which will be 
discussed in the next subsection. 

For future reference we note that the linear combination: 

Q ■= le~ 2p <p'g + IGe-^ILpF 

yields a new hamiltonian constraint in terms of the mass function: 

M = ± [e^{{GU p f - (0') 2 ) + J -^) (17) 

M. is called the mass function because it approaches a constant at spatial 
infinity, where it is equal to the ADM mass of the solution. In the presence 
of matter, Q can be integrated to yield a simple and physical constraint that 
relates the mass function at any spatial coordinate to the integrated matter 
stress energy: 

With appropriate boundary conditions, to be discussed in the next sec- 
tion, and corresponding Hamiltonian boundary term needed to make the 
variational principle well defined, the equations of motion for the momenta 
are: 

on, = -°"-(°py+^ e2p %+( NGu *y ( i9 ) 

Gtl p = (GNU p y-(a<t>y + ae 2 f^- (20) 

= (<Th(<j>wy + (Nu i ,y (2i) 
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4 Partial Gauge Fixing and Boundary Con- 
ditions 



As in [61 [8] we choose a partial gauge and boundary conditions that ensure 
PG metric at infinity and in vacuum. In [B] a timelike gauge fixing condition 
was imposed, and the Hamiltonian constraint was set strongly to zero. We 
choose instead to fix the spatial coordinate by adding the constraint: 

X = 0'-J'(0)~O (22) 

As usual, this is imposed weakly (after the Poisson brackets are calculated). 
In order for this to be a good gauge fixing condition, it must be second class 
with at least one of the two original constraints. In fact it is second class 
with the diffeomorphism constraint: 

{ X (A),^(r)} = J dxTj(<P) (V + ^Ma) (23) 

which is invertible for all field configurations on the space of test functions A 
that vanish on both boundaries. 

Consistency requires that x must be preserved in time, 

X = { X ,H} = (24) 

which gives us the condition on the lagrange multipliers: 

cxGUp = N<j>' (25) 

Since all the phase space variables besides II^ commute with x the Dirac 
brackets for p, Yl p , ip, IT^ are the same as the Poisson brackets. This means 
that we can set x an d T strongly to zero in the Hamiltonian, and obtain the 
correct equations of motion using the ordinary Poisson brackets. 

After some algebra one can show that the reduced Lagrangian is: 

C = J dx (U p p + n^) - J dx^-g - J dx(^M)' (26) 

where Q is defined in (1161) . We have added the boundary term needed for 
slicings that approach PG coordinates at both boundaries. The partially 
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reduced equations of motion are: 



V 



p = IG[^-) e - 2 m p + lGa-^- (27) 





5>+^) 8 + 2 ""^7 y ) (28) 

j(0) /i(0) 

= (so) 

In order to discuss boundary conditions we need to specify the functions 
j{4>) and h{4>) that essentially define the theory. In order to admit single 
horizon black hole solutions, j(<ft) must be monotonic and without loss of 
generality we assume it vanishes at = 0, which therefore corresponds to 
the singularity. For simplicity we assume j(4>) = <p a where < a < 2. It was 
shown in [T3J that this latter condition is necessary for the theory to describe 
black holes. We also assume that the coupling h(<p) takes the simple form 
h((p) = <f) b where b > 0. For spherically symmetric gravity in D = 2 + n 
dimensions, a = 1 — 1/n and 6=1. Under these circumstances suitable set 
of boundary and fall-off conditions are given as follows. At spatial infinity, 
we require: 



e 



-(700-0)' 



2p - m (i+oor 7 -)) 

GTl p -> v /2GM ooJ //(l + 0( 

O -> ^ (l + 0(</.-^ +a / 

V> -> o(^ fe/2 - eo °) 

-> o(0 6 / 2 - e -) (3i) 

where 7^ > 1 — a/2, and again and Coo are parameters. 

At the origin, = 0, the fall off conditions require a bit more care. It 
turns out that the slicings can be extended to the origin in the presence of a 
singularity by requiring: 

e 2p -> j(0)(l + O(0 W2+70 )) 
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^2GM j/l + O(<j> 1+ ^) 
O(0 2 - a+70 ) 

0^-1/2+6/2+6) (32) 

where 70 and e can take on any positive value. In the above, M is the mass 
of the singularity. Note that since <r(0) = 0, we are restricting the spatial 
slices to approach the same PG slice at the origin. This requires M to be 
independent of time, since the spatial slice does not evolve there. The slic- 
ings can, however, evolve arbitrarily at the horizon, so that these boundary 
conditions are suited to the study of horizon formation and evolution prior 
to singularity formation. 

It can be verified that these boundary conditions guarantee that there are 
no additional contributions to the surface terms in the action, the Hamilto- 
nian and Liouville terms in the action are finite, and that they are preserved 
by the dynamical evolution equations. In fact, it is the consistency of the 
latter that forces us to restrict to spatial slices that do not evolve at = 0. 

Ideally one would like to be able to find more general fall-off conditions 
that allow a description of not only horizon formation, but also singular- 
ity formation. That is, one should start with a spacetime and scalar field 
which are initially regular at r = 0. At late times, some time after horizon 
formation, one expects a spacelike singularity at r = 0. In principle, our for- 
malism should allow a description of such a transition by fixing the behavior 
of the constant of integration M (t) in Eq. (ll8p . which is in principle an ar- 
bitrary function of time. In our choice of parameterization, a non-zero value 
for M (t) signals the presence of the curvature singularity at the origin. At 
large times, as the scalar field falls to the origin, M (t) should approach the 
ADM mass. Unfortunately, we have been unable to find consistent boundary 
conditions that allow non-zero, time varying M (t). This is currently under 
investigation. 

5 A Canonical Transformation 

The simplicity of the partially reduced Hamiltonian is particularly striking 
after the following canonical phase space variables for the remaining gravi- 
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tational field: 

X(x) := e p 

P(x) := e- p GU p (33) 
The Poisson brackets now read: 

{X(x),P(y)} = G8(x,y) (34) 
The full Hamiltonian in this case takes a very simple and suggestive form: 

H(X, P, ^ H^) = [dx i-^rM' + oQ M + *l XP fJ** ) + / dx{^-M)' 

(35) 

where 

and Qm is defined in (111]) . 

Eq. (|35|) is the main result of the present paper. It has some remarkable 
features. First of all, the Hamiltonian has decomposed into three terms: the 
first term involves only the gravitational phase space variables, the second 
only the matter variables, while the third term provides the coupling between 
them. Moreover, the coupling term is very simple in that it is quartic in the 
dynamical variables. Note that this in part is a consequence of our partial 
gauge fixing: the function <p(x) is no longer dynamical but a fixed function 
of the spatial coordinates. 

The key operators that appear in the Hamiltonian are: K := X 2 /4, R : = 
P 2 - &g-,D := IF/2, (</>') 2 , (n^,) 2 , and U^'. Note that K, R and D obey 
an 5*0(2, 1) algebra, so that the operators D and K give rise to classically 
conserved charges under time evolution generated by R. This algebra is 
known to be anomalous at the quantum level and has been extensively studied 
using a variety of techniques [HI [T5] . 



6 Conclusions: Toward the Quantum Dynam- 
ics of Black Hole Formation 

Equation fl35l) provides an intriguing starting point for the study of the quan- 
tum dynamics of gravitational collapse, where one might follow a procedure 
similar to the one outlined in Husain and Winkler [B]: 
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1. Find a suitable quantum representation for H. 

2. Fix the lapse function thereby implicitly choosing a time coordinate. 

3. Find a suitable initial state |\&(to)) that solves the constraint equation 
and describes the quantum version of some (semiclassical) initial scalar 
field configuration. 

4. Evolve the initial state (numerically if the representation allows) via 
infinitesimal time steps via the equation: 

|tf(to + At)) = (l + iifAt)|tf(to)) (37) 

5. At each time step check whether the expectation value of the future 
null expansion operator vanishes somewhere on the spatial slice. That 
is one looks for a value of for which: 

(V(t)|M*AorM*)> = (38) 

where 9 + (x) is a suitably quantized version of the classical future null 
expansion [T6| [7J: 

9 + (x) = P-£ (39) 

One would expect that in different gauges x\ lOT would take on different 
values, but that the presence, or not, of a horizon would be invariant. 

The key step in the procedure for quantum evolution is to start with 
a suitable representation which allows the constraint to be solved and the 
Hamiltonian evolution to be evaluated with finite computational resources. 
In [7J the Hamiltonian constraint in the vacuum case was solved in the 
Bohr representation in terms of a set of eigenstates \M) of the operator 
At. The Bohr quantization scheme arises naturally in the context of loop 
quantum gravity (see [17] for a review) but has also been recently advocated 
as useful for singularity resolution in the quantization of purely geometrical 
variables [IS]. The states found in [7J had interesting properties, including the 
fact that they induced a specific (mass dependent) discretization of space ex- 
terior to the horizon. This spectrum approached the continuum far from the 
horizon as well as for macroscopic black holes near the horizon. It was argued 
that they provide a reasonable basis for the Hilbert space of physical states 
of an isolated black hole. A basis of mass eigenstates was also obtained in [S] 
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using a Schrodinger representation. In this representation as well interesting 
quantum properties were observed near the horizon. 

Irrespective of which representation is used, Eq. fl35|) in conjunction with 
the results of [TJ, [8] suggests a natural basis for studying the quantum dy- 
namics of black hole formation. The first step would be to find a complete 
set of eigenstates, which we will denote generically by \k), of the operator Q M 
which is the Hamiltonian density for the matter fields. Given the structure 
of Eq. (1351) one could then start with the following representation for the time 
dependent quantum states. 

\^(t)) = Y.C M ; k (t)\M;k) (40) 

M,k 

where \M, k) := \M) <g> \k) is a direct product of eigenstates of the two 
operators. Clearly the action of the first two terms in fl35l) is relatively 
simple to compute, and all the interesting dynamics will be generated by the 
interaction term, which is also relatively simple as a quantum operator. 

Although this procedure is highly non-trivial, it does on the surface ap- 
pear to be tractable. We are currently working on implementation via the 
Bohr representation. 
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